Interfacial Organization and Forces Arising from Epithelial-Cancerous Monolayer Interactions.

The interfacial interactions between epithelia and cancer cells have profound relevance for tumor development and metastasis. Through monolayer confrontation of MCF10A (nontumorigenic human breast epithelial cells) and MDA-MB-231 (human epithelial breast cancer cells) cells, we investigate the epithelial-cancerous interfacial interactions at the tissue level. We show that the monolayer interaction leads to competitive interfacial morphodynamics and drives an intricate spatial organization of MCF10A cells into multicellular finger-like structures, which further branch into multiple subfinger-like structures. These hierarchical interfacial structures penetrate the cancer monolayer and can spontaneously segregate or even envelop cancer cell clusters, consistent with our theoretical prediction. By tracking the substrate displacements via embedded fluorescent nanobeads and implementing nanomechanical modeling that combines atomic force microscopy and finite element simulations, we computed mechanical force patterns, including traction forces and monolayer stresses, caused by the monolayer interaction. It is found that the heterogeneous mechanical forces accumulated in the monolayers are able to squeeze cancer cells, leading to three-dimensional interfacial bulges or cell extrusion, initiating the p53 apoptosis signaling pathways of cancer cells. We reveal that intercellular E-cadherin and P-cadherin of epithelial cells differentially regulate the interfacial organization including migration speed, directionality, spatial correlation, F-actin alignment, and subcellular protrusions of MCF10A cells; whereas E-cadherin governs interfacial geometry that is relevant to force localization and cancer cell extrusion, P-cadherin maintains interfacial integrity that enables long-range force transmission. Our findings suggest that the collaborative molecular and mechanical behaviors are crucial for preventing epithelial tissues from undergoing tumor invasion.

Two-point optical manipulation reveals mechanosensitive remodeling of cell-cell contacts in vivo.

Biological tissues acquire reproducible shapes during development through dynamic cell behaviors. Most of these behaviors involve the remodeling of cell-cell contacts. During epithelial morphogenesis, contractile actomyosin networks remodel cell-cell contacts by shrinking and extending junctions between lateral cell surfaces. However, actomyosin networks not only generate mechanical stresses but also respond to them, confounding our understanding of how mechanical stresses remodel cell-cell contacts. Here, we develop a two-point optical manipulation method to impose different stress patterns on cell-cell contacts in the early epithelium of the Drosophila embryo. The technique allows us to produce junction extension and shrinkage through different push and pull manipulations at the edges of junctions. We use these observations to expand classical vertex-based models of tissue mechanics, incorporating negative and positive mechanosensitive feedback depending on the type of remodeling. In particular, we show that Myosin-II activity responds to junction strain rate and facilitates full junction shrinkage. Altogether our work provides insight into how stress produces efficient deformation of cell-cell contacts in vivo and identifies unanticipated mechanosensitive features of their remodeling.

Structure and Rheology in Vertex Models under Cell-Shape-Dependent Active Stresses.

Biological cells can actively tune their intracellular architecture according to their overall shape. Here we explore the rheological implication of such coupling in a minimal model of a dense cellular material where each cell exerts an active mechanical stress along its axis of elongation. Increasing the active stress amplitude leads to several transitions. An initially hexagonal crystal motif is first destabilized into a solid with anisotropic cells whose shear modulus eventually vanishes at a first critical activity. Increasing activity beyond this first critical value, we find a re-entrant transition to a regime with finite hexatic order and finite shear modulus, in which cells arrange according to a rhombile pattern with periodically arranged rosette structures. The shear modulus vanishes again at a third threshold beyond which spontaneous tissue flows and topological defects of the nematic cell shape field arise. Flow and stress fields around the defects agree with active nematic theory, with either contractile or extensile signs, as also observed in several epithelial tissue experiments.

Mechanical stress driven by rigidity sensing governs epithelial stability.

Epithelia act as a barrier against environmental stress and abrasion and in vivo they are continuously exposed to environments of various mechanical properties. The impact of this environment on epithelial integrity remains elusive. By culturing epithelial cells on 2D hydrogels, we observe a loss of epithelial monolayer integrity through spontaneous hole formation when grown on soft substrates. Substrate stiffness triggers an unanticipated mechanical switch of epithelial monolayers from tensile on soft to compressive on stiff substrates. Through active nematic modelling, we find that spontaneous half-integer defect formation underpinning large isotropic stress fluctuations initiate hole opening events. Our data show that monolayer rupture due to high tensile stress is promoted by the weakening of cell-cell junctions that could be induced by cell division events or local cellular stretching. Our results show that substrate stiffness provides feedback on monolayer mechanical state and that topological defects can trigger stochastic mechanical failure, with potential application towards a mechanistic understanding of compromised epithelial integrity during immune response and morphogenesis.

Drosophila mechanical nociceptors preferentially sense localized poking.

Mechanical nociception is an evolutionarily conserved sensory process required for the survival of living organisms. Previous studies have revealed much about the neural circuits and sensory molecules in mechanical nociception, but the cellular mechanisms adopted by nociceptors in force detection remain elusive. To address this issue, we study the mechanosensation of a fly larval nociceptor (class IV da neurons, c4da) using a customized mechanical device. We find that c4da are sensitive to mN-scale forces and make uniform responses to the forces applied at different dendritic regions. Moreover, c4da showed a greater sensitivity to localized forces, consistent with them being able to detect the poking of sharp objects, such as wasp ovipositor. Further analysis reveals that high morphological complexity, mechanosensitivity to lateral tension and possibly also active signal propagation in dendrites contribute to the sensory features of c4da. In particular, we discover that Piezo and Ppk1/Ppk26, two key mechanosensory molecules, make differential but additive contributions to the mechanosensitivity of c4da. In all, our results provide updates into understanding how c4da process mechanical signals at the cellular level and reveal the contributions of key molecules.

Being able to sense harm is essential for survival. Animals have to be able to tell the difference between a gentle touch and a dangerous pressure. They do this using nerve cells called mechanical nociceptors which switch on when the body feels a potentially painful pressure, such as a sharp object poking the skin. Once activated, the nerves send outputs to other parts of the central nervous system which coordinate the motions needed to escape the source of the pain. One popular model to understand harm-sensing is the larvae of fruit flies which automatically roll back and forth when they sense the pointy sting of a wasp. This process is initiated by sensory nerve cells called class IV dendritic arborization neurons (or c4da for short) which sit under the fly's skin. However, it is still not fully understood how these mechanical nociceptors detect the poking forces of the wasp's tail. To investigate, Liu, Wu et al. built a device that could poke sections of fly larvae under a microscope so they could see how different types of pressure affected the activity and shape of c4da cells. This revealed that c4da nerves were most sensitive to sharp objects that illicit a more localized force, which may explain why these cells are so good at responding to wasp attacks. Further analysis showed that this sensitivity was due to the high number of branches, or dendrites, protruding from the body of c4da nerves. Liu, Wu et al. discovered that the dendrites were coated in a touch-sensitive protein that can sense and amplify both squashing and pulling, resulting in a signal that activates c4da nerves to send outputs to other parts of the central nervous system. This mechanism increases the likelihood that a c4da cell will detect a mechanical pressure even if it is far away from the body of the nerve. These findings shed light on how sensory cells like c4da are optimized to carry out specific roles. This could be important for understanding other nerve systems which sense mechanical pressure, such as those involved in touch or auditory processes. However, further work is needed to see whether the molecules and mechanism identified by Liu, Wu et al. are also present in humans.

Spontaneous organization and phase separation of skyrmions in chiral active matter.

Skyrmions are topologically protected vortex-like excitations that hold promise for applications such as information processing and electron manipulation. Here we combine theoretical analysis and numerical simulations to show that skyrmions can spontaneously emerge in chiral active matter without external confinements or regulation. Strikingly, these activity-driven skyrmions can either self-organize into a periodic, stable square lattice consisting of half Néel skyrmions and antiskyrmions, where the in-plane flows display an antiferromagnetic vortex array, or undergo phase separation between skyrmions with different topological numbers. We identify that the emerging skyrmion dynamics stems from the competition between the chiral and polar coherence length scales dictated by the interplay of intrinsic chirality, polarity, and elasticity in the system. Our results reveal unanticipated topological excitations, self-organization, and phase separation in non-equilibrium systems and also suggest a potential way towards engineering complicated bespoke skyrmionic structures through manipulating active matter.

Active Nematic Flows over Curved Surfaces.

Cell monolayers are a central model system in the study of tissue biophysics. In vivo, epithelial tissues are curved on the scale of microns, and the curvature's role in the onset of spontaneous tissue flows is still not well understood. Here, we present a hydrodynamic theory for an apical-basal asymmetric active nematic gel on a curved strip. We show that surface curvature qualitatively changes monolayer motion compared with flat space: the resulting flows can be thresholdless, and the transition to motion may change from continuous to discontinuous. Surface curvature, friction, and active tractions are all shown to control the flow pattern selected, from simple shear to vortex chains.

Implementation of cellular bulk stresses in vertex models of biological tissues.

Vertex models describe biological tissues as tilings of polygons. In standard vertex models, the tissue dynamics result from a balance between isotropic stresses, which are associated with the bulk of the cells, and tensions associated with cell-cell interfaces. However, in this framework it is less obvious how to describe anisotropic stresses arising from the bulk of cells. In epithelia, such bulk anisotropic stresses could arise for instance through medial myosin fluctuations. Two recent publications-Tlili et al. (Proc Natl Acad Sci USA 116(51):25430-25439, 2019) and Comelles et al. (eLife 10:e57730, 2021)-have proposed different schemes to implement bulk anisotropic stresses in vertex models. Here we show that while both schemes transform in the same way under affine deformations, they lead to significantly different tissue dynamics. Our results are consistent with the interpretation that the Tilli et al. scheme describes bulk stresses that are uniform within each cell, while the Comelles et al. scheme corresponds to non-uniform bulk stresses. Finally, we wondered whether a standard vertex model can be fully expressed in terms of bulk cellular stresses alone. We find that, in general, neither scheme can mimic the vertex forces created by cell-cell interface tensions.

Universal Statistical Laws for the Velocities of Collective Migrating Cells.

Migratory dynamics of collective cells is central to the morphogenesis of biological tissues. The statistical distribution of cell velocities in 2D confluent monolayers is measured through large-scale and long-term experiments of various cell types lying on different substrates. A linear relation is discovered between the variability and the mean of cell speeds during the jamming process of confluent cell monolayers, suggesting time-invariant distribution profile of cell velocities. It is further found that the probability density function of cell velocities obeys the non-canonical q-Gaussian statistics, regardless of cell types and substrate stiffness. It is the Tsallis entropy, instead of the classical Boltzmann-Gibbs entropy, that dictates the universal statistical laws of collective cell migration. The universal statistical law stems from cell-cell interactions, as demonstrated by the wound healing experiments. This previously unappreciated finding provides a linkage between cell-level heterogeneity and tissue-level ensembles in embryonic development and tumor growth.

Pattern Formation and Defect Ordering in Active Chiral Nematics.

Many biological systems display intriguing chiral patterns and dynamics. Here, we present an active nematic theory accounting for individual spin to explore the collective handedness in chiral rod-shaped aggregations. We show that coordinated individual spin and motility can engender a vortex-array pattern with chirality and drive ordering of topological defects. During this chiral process, the stationary trefoil-like defects self-organize into a periodic, hexagon-dominated polygonal network, which segregates persistently rotating cometlike defects in pairs within each polygon, leading to a translation symmetry at the global scale while a broken reflection symmetry at the local scale. Such defect ordering agrees exactly with the Voronoi tiling of two-dimensional space and the emergence of the hexagonal symmetry is deciphered in analogy with topological charge neutralization. We calculate energy barriers to the topological transition of the defect ordering and explain the existing metastable states with nonhexagonal polygons. Our findings shed light on the chiral morphodynamics in life processes and also suggest a potential route towards tuning self-organization in active materials.

Collective dynamics of coherent motile cells on curved surfaces.

Cellular dynamic behaviors in organ morphogenesis and embryogenesis are affected by geometrical constraints. In this paper, we investigate how the surface topology and curvature of the underlying substrate tailor collective cell migration. An active vertex model is developed to explore the collective dynamics of coherent cells crawling on curved surfaces. We show that cells can self-organize into rich dynamic patterns including local swirling, global rotation, spiral crawling, serpentine crawling, and directed migration, depending on the interplay between cell-cell interactions and geometric constraints. Increasing substrate curvature results in higher cell-cell bending energy and thus tends to suppress local swirling and enhance density fluctuations. Substrate topology is revealed to regulate both the collective migration modes and density fluctuations of cell populations. In addition, upon increasing noise intensity, a Kosterlitz-Thouless-like ordering transition can emerge on both undevelopable and developable surfaces. This study paves the way to investigate various in vivo morphomechanics that involve surface curvature and topology.

Dynamic instability and migration modes of collective cells in channels.

Migrating cells constantly experience geometrical confinements in vivo, as exemplified by cancer invasion and embryo development. In this paper, we investigate how intrinsic cellular properties and extrinsic channel confinements jointly regulate the two-dimensional migratory dynamics of collective cells. We find that besides external confinement, active cell motility and cell crowdedness also shape the migration modes of collective cells. Furthermore, the effects of active cell motility, cell crowdedness and confinement size on collective cell migration can be integrated into a unified dimensionless parameter, defined as the cellular motility number (CMN), which mirrors the competition between active motile force and passive elastic restoring force of cells. A low CMN favours laminar-like cell flows, while a high CMN destabilizes cell motions, resulting in a series of mode transitions from a laminar phase to an ordered vortex chain, and further to a mesoscale turbulent phase. These findings not only explain recent experiments but also predict dynamic behaviours of cell collectives, such as the existence of an ordered vortex chain mode and the mode selection under non-straight confinements, which are experimentally testable across different epithelial cell lines.

Three-dimensional collective cell motions in an acinus-like lumen.

Collective cell migration is an essential process in embryo development, wound healing, inflammatory response, and cancer invasion. Although cell motions in two-dimensional (2D) monolayers have been studied previously, three-dimensional (3D) collective cell migration, which constantly occurs during embryogenesis such as the establishment of ducts and acini in vivo, remains elusive. In this paper, we develop a cell-based model incorporating cell mechanics and cell motility to address coherent cell motions in a spherical acinus-like lumen with different cell populations. It is found that the interplays between cell persistence, random fluctuation, and geometrical confinement may engender rich and novel migratory modes. In a 3D spherical lumen, two cells may undergo stripe-like or cross-circular coherent rotations, whereas multiple cells can form dynamic twisting or circulating bands, leaving sparse cells at the center or even a hollow cavity in the cell aggregate. The cell density is found to profoundly influence the collective cell migration modes. Our model can reproduce the fundamental features observed in experiments and highlight the role of mechanics in steering 3D collective cell dynamics during mammary acinar morphogenesis.

Dynamic Migration Modes of Collective Cells.

Collective cell migration occurs in a diversity of physiological processes such as wound healing, cancer metastasis, and embryonic morphogenesis. In the collective context, cohesive cells may move as a translational solid, swirl as a fluid, or even rotate like a disk, with scales ranging from several to dozens of cells. In this work, an active vertex model is presented to explore the regulatory roles of social interactions of neighboring cells and environmental confinements in collective cell migration in a confluent monolayer. It is found that the competition between two kinds of intercellular social interactions-local alignment and contact inhibition of locomotion-drives the cells to self-organize into various dynamic coherent structures with a spatial correlation scale. The interplay between this intrinsic length scale and the external confinement dictates the migration modes of collective cells confined in a finite space. We also show that the local alignment-contact inhibition of locomotion coordination can induce giant density fluctuations in a confluent cell monolayer without gaps, which triggers the spontaneous breaking of orientational symmetry and leads to phase separation.

A nonlinear poroelastic theory of solid tumors with glycosaminoglycan swelling.

Mechanics plays a crucial role in the growth, development, and therapeutics of tumors. In this paper, a nonlinear poroelastic theory is established to describe the mechanical behaviors of solid tumors. The free-swollen state of a tumor is chosen as the reference state, which enables us to avoid pursuing a dry and stress-free state that is hard to achieve for living tissues. Our results reveal that the compression resistance of a tumor is primarily attributed to glycosaminoglycan (GAG) swelling, and the compactness of cell aggregates is found to affect tumor consolidation. Over-expressed GAGs and dense cell aggregates can stiffen the tumor, a remodeling mechanism that makes the tumor with higher elastic modulus than its surrounding host tissues. Glycosaminoglycan chains also influence the transient mechanical response of the tumor by modulating the tissue permeability. The theoretical results show good agreement with relevant experimental observations. This study may not only deepen our understanding of tumorigenesis but also provide cues for developing novel anticancer strategies.

Activation and synchronization of the oscillatory morphodynamics in multicellular monolayer.

Oscillatory morphodynamics provides necessary mechanical cues for many multicellular processes. Owing to their collective nature, these processes require robustly coordinated dynamics of individual cells, which are often separated too distantly to communicate with each other through biomaterial transportation. Although it is known that the mechanical balance generally plays a significant role in the systems' morphologies, it remains elusive whether and how the mechanical components may contribute to the systems' collective morphodynamics. Here, we study the collective oscillations in the Drosophila amnioserosa tissue to elucidate the regulatory roles of the mechanical components. We identify that the tensile stress is the key activator that switches the collective oscillations on and off. This regulatory role is shown analytically using the Hopf bifurcation theory. We find that the physical properties of the tissue boundary are directly responsible for synchronizing the oscillatory intensity and polarity of all inner cells and for orchestrating the spatial oscillation patterns inthe tissue.

Collective dynamics of cancer cells confined in a confluent monolayer of normal cells.

Tumorigenesis often involves specific changes in cell motility and intercellular adhesion. Understanding the collective cancer cell behavior associated with these specific changes could facilitate the detection of malignant characteristics during tumor growth and invasion. In this study, a cellular vertex model is developed to investigate the collective dynamics of a disk-like aggregate of cancer cells confined in a confluent monolayer of normal cells. The effects of intercellular adhesion and cell motility on tumor progression are examined. It is found that the stresses in both the cancer cells and the normal cells increase with tumor growth, resulting in a crowded environment and enhanced cell apoptosis. The intercellular adhesion between cancer cells and normal cells is revealed to promote tumor growth and invasion. The tumor invasion dynamics hinges on the motility of cancer cells. The cancer cells could orchestrate into different collective migration modes, e.g., directional migration and rotational oscillations, dictated by the competition between cell persistence and local coordination. Phase diagrams are established to reveal the competitive mechanisms. This work highlights the role of mechanics in regulating tumor growth and invasion.